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I. INTRODUCTION 

Studies of crack propagation and stability have developed in two main 
directions. One, related to materials science, concerns studies of the 
hierarchy of microdefects, their nucleation, interaction and development in 
association with propagation of a main crack. The microdefects are very 
different by nature depending on material structure and method of observa- 
tion, specifically the level of magnification chosen. Using a progressively 
finer scale of observation, a hierarchy of defects can be visualized. For 
instance, the various elements of damage in the vicinity of a fatigue crack 
in ASI 301 stainless steel are shown in Fig. 1 [1], A zone of a large 
plastic deformation surrounding the crack tip appears under 25X (Fig. la). 
The randomly oriented, "turbulent" field of lines representing the locali- 
zation of deformation can be distinctly observed under 500X (Fig. lb). 
Elements of discontinuity constituting an essential part of overall defor- 
mation can be identified under 2O,O0OX (Fig. lc). Obviously, other details 
could be seen on intervening magnifications. Individual dislocations, the 
atomic structure etc can be observed under larger magnification. Which 
elements of this hierarchy of defects should be parametrized in order to be 
included into a quantitative model of a crack surrounded by damage? Ap- 
parently this question should not be addressed to the materials science 
only. The continuum mechanics, constituting the foundatiaon of the second 
main direction in studies of fracture propagation, is addressed as well. 

Conventionally in continuum mechanics a crack is considered as an ideal 
cut in an elastic, elasto-plastic or visco-elasto-plastic medium. The 
concept of crack-cut with associated surface energy was the first and a 
very important step in studies of brittle failure. It reflects some essen- 
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Morphology of the crack layer on different magnifications, 
a - General view of the fatigue crack layer at low magni- 
fication. 

b - The crack tip region. Extensive damage is seen around 
and in front of the crack tip. 
c - SEM picture of an element of damage from the area in b 
taken at 20,000x magnification. 
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tial features of fracture processes and has served as a solid foundation 
for many engineering and scientific applications. Following this approach, 
microdefects surrounding the crack are modeled as a plastic zone in a very 
general sense. Such modeling describes the macroscopic deformation and 
stress state reasonably well for plastic metals. It does not, however, 
describe the microstructure of plastic deformation within the "plastic 
zone" (see for example Fig. lb, c). For brittle materials such as ceramics 
and rocks the models of plasticity seem inadequate. 

Recent achievements in materials science challenge the continuum 
approach to model the fracture processes. Obviously the complexity of a 
hierarchy of interacting defects briefly mentioned above is the main 
obstacle . 

In order to simplify the picture we may examine a crack surrounded by 
damage under relatively low magnification (Figs. 2,3,4) [1-3]. These 
micrographs have been obtained from various materials: polystyrene (an 
amorphous polymer), polypropylene (a semicrystalline polymer) and stainless 
steel (a polycrystalline metal). The observations [1-3] indicate that 
under similar loading conditions the global geometry and the evolution of 
an array of microdefects surrounding a main crack have many similar 
features for various materials in spite of all the differences in molecular 
structure and morphology. 

Fracture propagation is usually an irreversable process. Hence the 
general framework of the thermodynamics of irreversable processes can be 
employed for modeling the phenomenon. 

A system of a crack and its surrounding damage is referred to as a 
crack layer (CL). The theory of crack layer propagation based on irreversi- 
ble thermodynamics has been proposed in [4,5, 6]. Supporting experimental 
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evidences were reported recently [1-3,7, 8 ]. 

This work gives a comprehensive presentation of crack-layer (CL) as a 
model encompassing essential features of fracture propagation in various 
materials. Relation of the CL model to applied fracture mechanics problems 
is discussed. 

II# Thermodynamics of an Elastic Medium With Damage 

Fracture processes consist of nucleation and growth of microdefects. 
For thermodynamic description of fracture one needs to introduce a list of 
parameters of state by incorporating a damage parameter. To my knowledge, 
L.M. Kachanov was the first one who did it explicitly [9]. Since then, 
numerous papers suggesting various damage parameters and constitutive 
equations for them have been published. Summaries of some of the proposed 
damage models can be found in [10,11]. Specific interpretation of a damage 
is vitally important for establishing a correspondence between experimental 
studies and a damage model. It is not so important for general thermody- 
namic analysis, which is presented below. However, it is always useful to 
have in mind a particular damage model. Therefore, in this section we 

introduce a damage parameter (P) following [12]. Surfaces supporting the 
discontinuities within an initially continuous solid are considered as the 
elements of damage (microcracks, crazes, shear bands, martensite transfor- 
mation). A damage parameter P is defined as a pairing of scalar damage 

density p (i.e., an area of discontinuity surfaces per unit volume [ P] * 

m 2 

— 3 ) and damage orientation parameter 0: P = (p, 0} . Thus, the following 

nr 

system of thermodynamic parameters of state is considered. 


{a, T, P) 


( 1 ) 



in 
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Here stress tensor o and the absolute temperature T constitute a conven- 
tional for an elastic medium list of parameters of state. This list is 
extended by adding the damage parameter P in order to describe the fracture 
processes. To deduce 'thermodynamic causes' of damage we derive the entropy 
production associated with damage. 

The local energy balance can be written as: 


u = 


= o : c - V • j 


( 2 ) 


Here u stands for the rate of internal energy density, c is a strain rate 
tensor, the product g:e represents the rate of work density, gives the 

rate of internal energy density due to heat transfer, ( stands for heat 
flux) . 

Considering small deformations, we decompose the total strain tensor 
into perfectly elastic e (thermodynamically reversible) and nonelastic 
(irreversible) parts, i.e., 



(3) 


The work done on nonelastic deformation is spent partially on damage 
and partially by being converted into heat. It can be expressed as 
follows: a a:e^ is a part of the irreversible work associated with 
damage nucleation and growth and (l-a)a:g^^ is converted into heat, (a is 
a phenomenological coefficient). 

Let us consider the left part of equation (2). Conventionally, the in- 
ternal energy density u consists of Helmholtz free energy density f and the 
entropic part Ts: 
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u * f + T s 


(4) 


consequently: 


u - f + Ts + ts (5) 

According to the basic concept of irreversible thermodynamics, the time 
rate of changes of the entropy density s can be decomposed into two terms: 

s = Si + s e , (6) 

where stands for the entropy production due to irreversible processes, 
and s^ is the entropy density rate due to exchanges with the surrounding by 
heat and other kinds of energy. In the equilibrial thermodynamics the 
entropy Increment As is defined as the ratio: heat/ temperature . In non- 
equilibrial thermodynamics the heat exchange with the surrounding is usual- 
ly assumed to be equilibrial # Using this assumption we Introduce the 
equilibrial entropy rate s^ in the following form: 

* - - S-Cjjr j Q ) + ^-2 o:c (i) + Asp (7) 

e i ~ T — ~ 

Here the first term represents the exchange entropy rate due to the entropy 
s 1 0 

flux j * j , the second term reflects the entropy increases due to the 
heat generated by the irreversible work 0:e ^ and the third term reflects 
the rate of the entropy changes due to the localized transformation of the 
termodynamic state such as cracking, crazing, shear banding, etc. As(g,T,o) 
stands for the difference between the entropies of damaged and undamaged 
matter. 

Since the stress tensor and the absolute temperature constitute the 
conventional part of the list of parameters of state, it is convenient to 
use Gibbs free energy density g. We define g as the difference between 
Helmholtz free energy density and the density of the work done on elastic 
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deformation : 


g 


- „• < e > 


O • £ 


( 8 ) 


Obviously, 


f - oic e - g + o:c^ e ^ (9) 

Substituting (3)-(7) into (2), solving energy balance equation (2) with 
respect to the entropy production s^ and taking into account (9) one an 
find: 

Ts i = a a:e (l) - g - 2: ? (e) - sT - ASp - y jV (io) 

For further transformations, we decompose the rate g into two terms: 

g(o,T,P) = Ag(a,T,0)p + Tr(o,T,P) (11) 

where Ag is the difference between the Gibbs free energy densities of 
damaged and undamaged matter and it is the elastic potential energy density. 

The assumption of local equilibrium yields the following conventional 
consitutive equations: 


3tt 

St |q 
p 


s 


( 12 ) 


3tt _ (e) 

9o|p 

T 


(13) 


9 


We also make use of conventional relationship between the densities of 
the enthalpy h, Gibbs free energy g and entropy s in the following way: 


Ag + TAs = Ah(g,T,0) 


(14) 


Then substituting (11)-(14) into (10) and having defined 1 P =Ahp the 


following expression for the entropy production can be obtained: 

Ts. =ao:^^ - +T1 - ) u _ Jl 


dP 


p - Y j • VT 


( 15 ) 


The last term in (15) describes the entropy production due to heat 
transfer which is out of the scope of this paper. In order to concentrate 
attention on the damage process (P), we assume an isothermal condition and 
homogeneity of temperature field (V T = 0). Under these conditions the 
entropy production (15) reduces to two terms, the first ao : is in- 

directly associated with damage and therefore is nonzero only within a 
region where P^O. The second term represents the entropy production due to 
damage growth (P) directly. Therefore, 

Ts. - a a:e (i) - ^ — P (16) 

i ~ ~ 3P 

The rate of damage growth P may be considered as a thermodynamic flux. 
_ 9 (h + tt) 

Then ^ is the reciprocal force. The physical interpretation of 

this force can be done using a simple example. If one visualizes the 
damage as a field of microcracks with microcrack density p then the rate of 
the potential energy density ti(o,T,P) with respect to p is always negative 
[12], i.e. 

9tt >0 
9p|g 
T 

It means that the contribution of the potential energy density change into 
the thermodynamic forces is always positive. It can be shown for the same 


1 See also equations (18) and (44). 
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conditions that the enthalpy increment Ah as well as the entropy increment 
As due to cracking is always positive, i.e., Ah >_ 0 [12]. Thus, the thermo- 
dynamic force - ^ reciprocal to the rate of damage P results from 

dr 

the competition between the driving (- ~ > 0) and the resisting (- ^ < 0) 
parts. 

To emphasize the importance of the enthalpy increment for the thermody- 
namics of failure we introduce the specific enthalphy of damage Y with 

2 

respect to the unit of damage chosen (p ^s): 

m - 5 

Y*(o,T,0) d = f Ah(Jtif 2 ) (18) 

The establishment of a constitutive equation for damage growth P would 
be the most desirable goal. This requires an effort to specify the damage 
parameter P and a wide program of microscopical observation of damage 
evolution under various loading conditions. However, at the present time 
we do not have sufficient experimental data to construct a constitutive 
equation for damage growth on a sound foundation. An attempt to do it, 
although useful, would be just a speculation with many adjustable para- 
meters Involved. In addition, for crack propagation studies, integral 
characteristics of damage surrounding the crack rather than details of the 
damage distribution, are Important. Therefore, following the spirit of 
thermodynamics we introduce average characteristics of the entire damage 
zone Including the crack (CRACK LAYER) and constitutive equations for the 
crack layer without reference to a constitutive law for local damage (P) . 
Obviously, the constitutive equations for the crack layer could be deduced 
from a constitutive law of local damage. The absence of this law prohibits 
the stated deduction. The crack layer approach presents a valuable 
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solution to this problem. 

.III. The Concept of CRACK LAYER 

We consider a system of a crack and the surrounding array of micro- 
defects as one macroscopic entity, crack layer (CL). As an illustrative 
example the trace of fatigue CL propagation in polystyrene (2] is shown in 
Fig. 5. 
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Fig- 5 

The vertical markers indicate the positions of the crack tip corresponding 
to various numbers of cycles. The damaged zone expands in a self similar 
fashion. The observation of damage (crazes) in the vicinity of the crack 
tip [2] suggests that the damage distribution appears as a manifestation of 
actual stress field. Consequently, the similarity criteria for damage 
distribution are the same as that for the stress field. 

Formally, a crack layer is described as zone within which the damage 
density p is positive, i.e., 

V L = {(x^ x 2 ):p(x 1> x 2 ) > 0 } (19) 

It is worth noting that a certain level of damage could exist independently 

i 

of crack propagation. In such a case, a level of reference damage density 

p should be determined. The CL is then defined as a zone V T within which 

o t. 


Number of Cycles 
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To describe CL propagation, we distinguish active and inert zones (see 


Zone hr 


stress concen 


tration. The zone V 


consequently the process of damage growth practically stops 


Thus, the inert zone appears as a trace of the active zone propagation 


The active zone boundary ?V can be represented as consisting of the lead 


of the active zone are shown 
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IV. Kinematics of CL Propagation 

The mechanism of crack layer propagation can generally be described as 
follows. At first, damage nucleates and grows to a critical level at which 
the local instability conditions are met [12]. Then a crack appears within 
the damage zone. The crack creates a stress concentration which 
intensifies the processes of damage growth. The damage density in front of 
the crack reaches a critical level which leads to a crack extension and so 
on . 

The critical level of damage relates to the stresses through the condi- 
tions of instability [12]. Since the stress distribution in the vicinity 
of a crack has an invarient shape with respect to crack length (the crack 
length just scales the stresses by a stress intensity factor K) , the criti- 
cal level of damage is maintained constant during the crack propagation. 

The crack layer propagation can be visualized as an active zone (damage 
distribution) movement. The latter can be decomposed into translation and 
rotation as a rigid body and an active zone deformation (i.e. damage dis- 
semination). Considering homogeneous deformation only, one can express the 
rate of the damage parameter £ resulting from the active zone movement 
without explicit time dependency in the following form: 


P = v 


r tr . 

6 P + ui 


o P + e 


«“ P P + d: 6 dev p 


( 21 ) 


Here v, u), e, <! stand for the rates of translation, rotation, isotropic 

expansion and deviatoric deformation, correspondingly, and 6 tr , 6 t0t , 6 eX * ) , 
dev 

6 are the operators of translational, rotational, expansional and devia- 


toric transformations 
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V. Global Entropy Production Due to Cr ack Layer Propagation 

The global entropy production is the integral over the entire volume of 
a solid V from the entropy production: 


•global 

b i 


p 

s. dV 

J 1 

V 


( 22 ) 


Since the entropy production due to damage growth (16) is nonzero only 
within the active zone, the integral in (22) is reduced into the integral 
over V^. Substituting (21) into (16) and integrating we obtain, 


Tg global = ^ 4 v-X tr 


+ urX 


rot 


+ e-X 


exp 


+ d:X 


dev 


(23) 


Here 0 is the part of irreversible work within the active zone spent on 
damage , i .e . , 


D = a [ o : dV = + 6 


(24) 


where * f 0:g^^ dV is the rate of the total work dissipated within 

Jv A 

and the rate of heat radiation Q *J j^.ndA. The thermodynamic forces 
translational X tr , rotational X rot , expentional X ex ^ and deviatoric X^ ev 
can be presented as follows: 


•tr 


f 9 < h + 

k 9P 


^ 6 tr P dV 


(25) 


x rot = f 90L+JI ) 6 rot p dy 


V 

A 


3P 


(26) 


exp 


_ f d (h + 7 

" Jv 3P 


jr) exp 


6 P dV 


(27) 
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V dev _ f 9(h + in ) x dev„ 

X * J gp § P dV (28 ) 

V A 

VI. Thermodynamic Forces for the CL 

In order to calculate the forces we consider the enthalpic (resisting) 
and potential energy release (driving) parts separately. Let us start with 
the driving parts. 

1 . Calculation of the potential energy release rates . The following 
transformation can be applied to the elastic potential energy density 
71 ( 0 ,T,P) using arbitrary operator n 6" and assuming homogeneous temperature 
field. 


9tt r 9 tt - 

_ 6P , - 53:6s 


(29) 


Therefore 


| |j| 3p dV - | [6. - g:6o]dV <30) 

V A V A 

Now one can use (30) for various operators. 

(a) The operator of an infinltisimal translation in k-th direction 
(opposit to the active zone advance) is 


tr _ _9_ def ^ 

K " 9 x r " K 


Substituting (31) into (30) we obtain 



(3 k " - 35 ; V ,dv 


(31) 


^ 6 tr P dV 
9P K 


J 


t 


(32) 
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Using the constitutive equation (13) ( _ ft? . = _ C ii), expression (8) written 

9o . . J 

e ij e 

in the following form Tr * f (e ,T,P) - 0 (f is the elastic strain 

energy), equations of equilibrium in the absence of mass forces: 


3 o - 0 
j ij 


and the definition of a small deformation 


f = X (u. . + u. .) 

IJ 2 1,J J,1 


the expression within the brackets in (32) can be simplified: 


« So ij * [ a U ij 1 ij Vi; 


V )* ' °i_j u 


ij i,j< 




Substituting (35) into (32) and using Green's theorem we find 


Ip 6 r T p dv “ J < * K = 1,2 


Where stands for the conventionally used energy release rate with re- 
spect to crack extension in the tangent direction: 


J 1 " P "l ‘ °ij n j U i’l ldT 

3V. 
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If one chooses a closed contour F outside of the active zone V^, then the 


integral (37) is identically equal to zero since the left part of equation 

(36) vanishes (see Fig. 7). The path independency of (37) is direct 

consequence of this statement. Therefore, instead of 3V^ ■ ^ 

in ( 37 ) , one can choose a part- 

ially arbitrary contour r * 

r (t) ur* (Fig. 7). The integral 

in (37) is invariant with re- 
ft 

spect to an arbitrary path F • 



(38) 


Fig. 7 

b) • The operator of an infinitesimal rotation is 

( correction for tensor components \ 
associated with rotation of coordinate I 
system spinorlal terms. / 

where is the alternating symbol with components equal to + 1 or - 1 

for even and odd permutations of 1, 2, 3 respectively and zero for all 
other combinations. 

Substituting (38) into (30) and providing transformations similar to 
that in (35), one can obtain: 


9 IT rTOt 

3p P4V 


m 


( 39 ) 


V, 


Here L is a pseudovector (m * 3 for 2-D problems) 
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L 

m 


f 


K£m 


! lY/ - 


3V. 


. . n . u . v 

ij J i.K 


<ViVK> d 


(40) 


For an isotropic medium the path independency of L remains similar to that 

m 

of for homogeneous medium* Consequently, the partially arbitrary con- 
(t) * 

tour I’* F Ul’ (see Fig. 7) can be used instead of ?V in (40). 

A 

c). The operator of an infinitesimal isotropic expansion is 


<S CX P = - 


K K 


(41) 


Substituting (41) into (30) one can obtain, 


6 exp PdV * M 


where M- f l e * °i j u t f « 3 dr 

3V* 


(42) 


(43) 


The M-integral (43) possesses the path-invariancy for linear medium only. 

The and M integrals represent the potential energy release rates 

with respect to the translation, rotation and isotropic expansion of the 

active zone. In a similar fashion, a second rank tensor N can be intro- 

KS. 

duced as the potential energy release due to deviatoric deformation of the 
active zone 16]. N ^ is particularly sensitive to deviations in the 
stress and strain fields near the crack tip. In contrast to J, L and M, N 
does not possess path independency. (See some historical remarks about J, L 
and M in the appendix.) 

Notably, the path independency of J, L and M renders their evaluation 
more convenient, it has no physical significance. For instance, for non- 
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homogeneous, anisotropic and nonlinear medium the equations (37), (AO) and 
(43) express the same potential energy rates as discussed above, but the 
integrals are not path independent. 

2. Calculation of the resistance parts of the thermodynamic forces. In 
order to calculate the enthalpic part of the thermodynamic forces one needs 
to specify the damage parameter. We make use of the damage parameter dis- 
cussed in Section II, i.e. P - { p , 0}, where p is the scalar damage 
density and 0 stands for an average damage orientation. If the orientation 
of an element of the damage (a plane microcrack, for instance) is described 
by a vector attached to the center of the element, then the translation and 
isotropic expansion transform the damage density P only and do not affect 
the orientation 0 . The rotation and the deviatoric deformation affect 
both the damage density and the orientation. 

The change of the enthalpy h due to damage according to (7), (11), (14) 
and (18) is given by: 

6P = Y*(o,T,CI)6p + p|j 60 (44) 


Here Y is the specific enthalpy of damage introduced in Section II. The 
general expression (44) can be used to derive the resistance parts of the 
thermodynamic forces. 

(a) Substituting (44), the expression (31) for the operator of trans- 
lations into the first term of (25) and keeping in mind that a translation 
affects the damage density only, one can find: 


f 3h 

J V 9P 

A 


6“P dV 


* 

Y 3„P dV = 

h K 

A 


* 

Y 


- 

j 


Pn K dr y \ 


3V 


( 45 ) 
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where R 


K 


P n^dT is defined as the translational resistance moment. 


3V, 


Let us consider a smooth crack trajectory 
(Fig. 8) and a local cartesian coordinate 
system with the origin at the crack tip 
and the axes ox^ directed along the tan- 
gent to the crack trajectory. Then, the 
crack speed vector v has only one nonzero 
component v ^( v 2 ~ 0)- Accordingly, only one component of the translational 
resistance moment is of interest. That is 



R i = - 


pn^dT = - 


3V. 


' I dr + 

r 1 

,(±) 


Pn^dF - 

w 


pn 1 dF = 
.(t) 


(46) 


pn l )<ir = l 


(p - P Q )dx 2 

.(t) 


where n v ~ y stands for a unit vector normal to 3V which has the positive or 
negative projection on the unit tangent x((n.t) ^ 0), respectively (Fig. 9).. 



Fig. 9 

and regular p r (x 2 ) densities, i.e., p(x 2 ) 

Q 

ly. R" is defined as the core of damage. 


The path of integration in (46) 
(the trailing edge) contains a 
singular point (the crack tip) 
where the damage density may be 
singular. Thus, we describe the 
damage density p(x 2 ) along 
as the sum of a singular R^ 6(xj) 
- R^ 6(* 2 ) + p T ( x 2 ^ » respective- 
substituting p(x 2 ) into (46) and 


integrating we find: 
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\ - (Rj + Rj] (46a) 

(b) In a similar fashion substituting (44) and the expression (41) for 
the isotropic expansion operator into the first term of (25) we obtain. 


|S «“ P P dV - - Y* 


X A° dV ' - YR „ 


V. 


(47) 


Here the expansion resistance moment R q is defined as: 


R = - 
o 


*kV' dv = - 


t» K (x K P) “ Pi dV = 


3V. 


X K pn K dF + j ° dV 


(48) 


In a case when the reference damage density P q » 0 the first term on 
the right hand side of (48) can be neglected. Indeed, 


j 

3V. 


dr = 


x K P o h K dr + 


X K Ph K dr ~ 0 


,(*) 


.(t) 


Since x^~0 and n^ - 0 along the trailing edge Therefore, 


R = p dV = <p> *A 
o 


(48a) 


where <p> stands for an average over damage density and A is the area 
of V A . 

(c) Substituting (44) and the expression (38) for the operator of 
rotation into the first term of (26) we obtain 


-I 


I C ’ 


- y*R - U R° 
m dO m 


( 49 ) 


Where R and R stand for rotational resistance moments, 
m m 
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Assuming that the enthalpic barrier does not depend on the orientation 

9-y* 

of damage (a reasonable assumption for an isotropic medium) i.e., — = 0 , 

aO 

then (49) can be reduced into 

- f S 6 r0t P dV - - y*R ( 50 ) 

J 9P m m 

V A 

For deviatoric deformation of the active zone a similar expression can 
be obtained: 


J 

V. 



P dV 


* 

- y 


(51) 


where 


IV A 


(1 - 




<V* 


x t 3 t )pdV 


( 52 ) 


and stands for Kroneker's symbol. 

Finally, summarizing the results of this section one can express the CL 
thermodynamic forces in an index form, which is convenient for practical 
application. 


*» - (J r + y'Rj) , 

(53) 

X™* - a 3 + W ot ), 

(54) 

X exp - (M + Y*R ) » 

O 

(55) 

Xjj ev - (N k - Y\ £ >. - 1,2 

(56) 


It should be noted that the singular damage density R^ 6 (x 2 ) contributes for 
translational resistance only, it produces no effects for R q , R ^ ot or 

\r 

VII. Rectilinear Crack Layer Propagation. Single Parameter Model 


To analyze the constitutive laws of crack layer propagation we start 
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from the simplest case when all degrees of freedom are frozen except one. 
We assume that crack layer propagation appears as a translation of the 
active zone along a rectilinear path with neither deformation nor rotation. 
Therefore, * 0; e = 0; ^ » 0 and * 0 (in terms of the coordinates of 

L 

Fig. 8). Since the crack trajectory is a rectilinear v^ » & where £ is the 
crack length. Using (53) the global entropy production (23) can be rewrit- 
ten as follows: 


T gglobal 


D + 



(57) 


According to the second law of thermodynamics, the entropy production 
(global as well as local) is nonnegative and equals zero for reversible 
processes. This does not mean that the processes causing negative entropy 
production not take place at all. Such processes are well known in 
chemical thermodynamics, [13,14]. The negative entropy producing processes 
may occur if some other dissipative processes produce a sufficient amount 
of entropy to make the total entropy production nonnegative. In this case 
the rates of entropy consuming processes are controlled by other sources of 
entropy production. The constitutive law can therefore be obtained from 
first principles. This is the case for the crack layer propagation. 

In order to derive the law of rectilinear CL propagation, we need to 
analyze the stability of the CL. Assuming that there are no other sources 
of dissipation except of CL growth, i.e., D = 0, the global entropy produc- 
tion takes the form 

Tg .global = ^ - y\) (58) 


In classical thermodynamics, the instability condition is the preroga- 
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tive of the second law. There are no disagreements about the criteria of 
•instability of an equilibrial state. However, it is not so clear what in- 
stability criteria should be accepted for an irreversible process. In this 
study we make use of the "universal criterion of evolution" which has been 
recently proposed and successfully applied to various irreversible process- 
es [14]. Since the entropy production is a bilinear form of thermodynamic 


fluxes j and reciprocal forces X^: 

s i- ( y k -\ • 

the time rate of the entropy production can be naturally decomposed into 


two terms: 


d§. d.S. d S. 
1 = J i , x i 

dt dt dt 


(60) 


where 


and 


fifi . ,, 

dt < k A dt 


d S. dX. 

4 _ k 

dt ( k ) k 


(61) 

(62) 


The criterion of evolution [14] states that the following inequality 


d S 

x l „ _ 
— 1 ° 


always holds true for stable processes. The equality is met for either a 

stationary process or a critical situation when a sudden (uncontrolled) 

2 . 

transition becomes possible. Analysis of the second variation 6 s^ is 
necessary to distinguish these cases. 

Let us apply this criterion to the entropy production due to crack 


layer extension (58): 
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. • global 

T -ar ■ * 2 • It (j i * * ,R i> 



< 0 


( 64 ) 


Indeed, the crack length l Is the only variable, since o is constant. is 
constant according to our assumption (there is no deformation of the active 
zone) . 

dj 

From (64) one can conclude that -jg- should be < 0 for a stable crack 
growth . 

Two types of crack configuration can be distinguished from the stabili- 
ty view point: stable (^ < 0) and unstable (jjj >0). Illustrative ex- 
amples of these two distinct types are shown in Fig. 10. 



Crack layer exhibits different behavior for stable and unstable configura- 
tions . 

1 ) CL propagation in a stable configuration 

For the configuration shown in Fig. 10a the energy release rate J is 

p2 J T 

given by Consequently ~ is always negative, i.e., the configuration 

is always stable. Then, the requirements of the second law (58) is the 
only "controller" of CL propagation. For J > y Rj the requirement (58) is 
met and a crack layer is "allowed" to grow with undefined speed i > 0. To 
specify the speed, we consider stationary CL growth. The condition of 
stationarity , i.e., the equality in (63) yields 

J 1 " 


(65) 
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Then, the applied load F controls the crack propagation velocity. Indeed, 

* 

from (65) and constancy of y one can deduce: 


d 

3t 


(J. 


y*Rj 


— F ♦ — £ 
3F h 3£ 16 


( 66 ) 


Substituting J 


into (66) we obtain 


* 

1 « 



(67) 


The condition (65) is obviously in agreement with the principle of 
minimum entropy production [21]. 

If J < Y the CL growth (i.e. £ > 0), consumes entropy (negative 

entropy production). According to the second law, it is possible if the 

• * 

dissipation D is sufficient to compensate the negative term £(J^- y R^). 
Applying the principle of minimum entropy production (the minimal value of 
^global zero), we find: 

f> ♦ i(J - y*R x ) - 0 (68) 

from which 


l 





Summarizing the cases discussed, one can write a general expression for 
a stationary CL propagation in a stable configuration: 


£ 


f'Rj* Jj if J i * Y * R i 

SJ 

- J i ■ ^* R i 

3£ 


Dissipation controlled 
process. 

(69) 

Equilibrial crack layer 
growth. 
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2 ) CL propagation in an unstable configuratio n 

For the configuration shown in Fig. 10b the energy release rate J * 

0 7T^' pvj O^TT 

— — . Consequently ^ * -g— is always positive, i.e. the configuration is 
always unstable. The limitation introduced by the second law prevents the 
crack from avalanche-like propagation, i.e., the crack cannot grow if 
TS« l0bal - i(Jj - Y*V < 0 

which yields (for l > 0) the following requirements for stable CL growth: 

J l - Y*R < 0 (70) 

Therefore, a slow CL propagation, which is compatible with (70) can occur 
due to the dissipative term t) only. Therefore, if one accepts the principle 
of minimum entropy production, then equation (68) holds true for an un- 
stable configuration as well as for a stable one. 

* 

When approaches y the requirement of the second law is met. Then 

for an unstable configuration the crack propagates avalanche-like if * 
* 

Y V 

Summarizing the results, one can write: 


* • »!>, if J 1 1 ** R 1 


Undefined if 


Dissipation controlled 
process. 

Transition to the (71) 
dynamic process. 


Although both characteristics 0 and y can be studied using either stable 
or unstable configuration, the first is more convenient for analysis of the 
dissipation 6 while unstable configuration is preferable for evaluation of 
y*R^. The later can be obtained from the conventional fracture toughness 
test and additional microscopic studies of damage distribution along 
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The rate of dissipation within the active zone t> can be experimentally 
measured as the difference between the rate of total dissipated work 0^ 
and the heat Q radiated (see 24). 

3) Example 

Slow CL propagation in unstable configuration (equation 71) is consi- 
dered as an illustrative example. Since, we do not have experimental data 
on heat radiation, we assume that the rate of dissipation 6 is proportional 
to the dissipated work The coefficient of proportionality apparent- 
ly depends on mechanisms of dissipation. One would expect the coefficient 
to be dependent upon strain rate, temperature and a characteristic time of 
the process. An evaluation of the irreversible work due to an array 
of crazes constituting the CL active zone has been done in [8]. The evalua- 
tion is based on a new CL stress analysis [15] and an experimental CL 
characterization [2]. These results suggests that is proportional to 
the product J<d>, where <d> stands for a characteristic size of V^. Thus, 
the rate of dissipation 0 can be written as follows: 


6 = B<d>J (72) 

Here g is a phenomenological coefficient with dimension [B ] = sec *. At 
fixed temperature B depends on the strain rate. The latter can be 
expressed in terms of the rate of applied load and a dimensionless crack 
propagation rate £/<d> [8]. 

Substituting (72) into (71) we obtain 


l - 


J if j<y*r, 
y*r x - j 1 

Undefined if J “ Y*R^ 


(73) 


* 

The relationship (73) is schematically represented in Fig. 11. One can 
distinguish three stages of slow crack propagation in unstable configura- 
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tion. The first stage corresponds to 



Fig. 11 


crack propagation through previously 
damaged material. A new damage 
is not being developed at this 
stage. It means that the re- 
sistance moment mainly con- 

£ 

sists of the core of damage R^ 
which is small in comparison 

with a developed CL resistance 
R^ . Therefore, the rate of 
crack acceleration with respect 
to the energy release rate 
is relatively high. This is in- 
dicated by the slope of the 
stage 1 portion of the curve 


in Fig. 11. 

The second, intermediate, state is characterized by monotonic growth of 

damage accompanied by crack propagation. It is reflected in a monotonic 

(approximately linear) increase of the translational resistance moment R^ 

and the characteristic length <d> with increasing of J^. In this case the 

translational driving force X tr is maintained approximately constant and 

the rate of crack propagation (73) can be approximated by a power type 
2 4 

equation l ~ J (i.e., K). 

The third, subcritical stage of CL propagation is characterized by 
deceleration of the resistance moment R^ with respect to J. This results 
in the translational force « '/‘R^ - J approaching zero, which corres- 
ponds to the critical state (the asymptote in Fig. 11). When y » 0 

the requirement of the second law (58) is met and instability becomes per- 
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missible. Thus, at the end of the third stage a slow crack growth trans- 
forms into uncontrolled (avalanche-like) crack propagation. 

Obviously, the value of J^at which the translational force Y*R^ - 

equals zero corresponds to the critical energy release rate J (or G ) in 

c c 

conventional fracture mechanics, i.e., 

J - Y* R lc (74) 

lc 

is a parameter which can be experimentally evaluated using fracture 
* 

mechanics, y and can be measured by materials science methods. Thus, 
(74) suggests a link between the two approaches. It is discussed in more 
details in section IX. 

The simple model described above generally predicts the shape of J vs. 
curve. However, it does not describe crack deceleration phenomena, 
history dependency of J^, etc. 

VIII. Rectilinear CL Propagation. Two Parameter Model 

Limitations of the previous model appear due to the employment of a 
single parameter only (the crack length ). A natural way to overcome the 
limitations is to use an additional degree of freedom offered by the CL 
model. Following this idea, we consider crack layer propagation by trans- 
lation along the rectilinear path and isotropic expansion of the active 

zone. Similar to the previous case, one can write: v^ - £ ; v^ * 0 (recti- 

A i 

linearity); u) ■ 0 (no rotation), e » 1/2 (— + ~) (isotropic expansion 

W a 

from the crack tip as an origin), d * 0 (no deviatoric deformation). 

The global entropy production (23) now can be rewritten as follows: 


T g global 


l (Jj - Y*Rj) ♦ e (M - y*R 0 ) (75) 


D + 
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The criterion of stability (63) for the crack layer only (ignoring D in 
<75) as we did in the previous example) can be written as follows: 


d s? lobal 

x 1 

dt 


co 2 . i; ■livr.'y . 


3(M- Y *R„) 


3£ 


] ♦ (e) 


3 £ 

2 3(M-y*R o ) 
3e 


3e 


< o 


Since £ and e are independent variables the above expression yields: 

< 0 
3 £ 

3 (M-y*R ) 

O 

< o 

3e 


3(J r Y*R 3 ) 

3£ 


*(J-y*R,) 

H - 


3e 


3 (J-Y‘R-1 ) 3(M-y*R 0 ) 

H- - ♦ - 

Be 3 1 


) 


5(M-y*R 0 ) 

3 £ 

3(M-Y«R 0 ) 

3c 


( 76 ) 
< 0 


Using the same argument as in the single parameter model, one can conclude 

* * 

that both the expansional (M - y R q ) and translational (J^ - yRj) forces 
are always nonpositive for a stable crack layer propagation. Therefore, 
two types of crack layer configurations (stable and unstable) can be dis- 
tinguished. Incorporating the of principle of minimum entropy production 
one can analyze stationary crack layer propagation following the formalism 
used in the single parameter model. 

The stationary crack layer propagation controlled by dissipation is 
described by the same equation (68) for both stable and unstable configura- 
tion (compare (69) and (71). For this reason, below we consider the dis- 
sipation controlled propagation only. 
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The principle of minimum entropy production for slow CL propagation can 
be expressed in the form: 

TSflobal - 6 - l (y\ ~ - e(y*R o -M) - 0 (77) 

It implies, that the dissipation 0 is distributed between two entropy sinks 
associated with two independent degrees of freedom: £ and e. As in the 

previous case we assumed the rate of dissipation 0 to be expressed by (72). 
In this case two parameters 3^ and S 2 substitute 3 to describe the distri- 
bution of the dissipation between two degrees of freedom. Then (77) yields 


. B <d >J 

a - 4 = — 


3 2 <d > J 

Y*R - M 
o 


(78) 


(79) 


Since the isotropic expansion is the only source for the resistance 
moment changes (it affects only the regular part R^) one can write 

Rj - R*(0)e (80) 

Accordingly, equation (79) can be converted into an equation for R* evolu- 
tion using (80): 


K 


B 2 R^(0)<d>J 1 

y\ ' M 


(81) 


To analyze the changes of r£ with crack length one can calculate 
by taking the ratio of (81) and (78): 


dRi 

d£ 



r Y * R 1 - J 
k R x y *R “ M 


( 82 ) 
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where k > j, is a phenomenological coefficient. It can be shown that J 
* * 

approaches y R_ faster than M approaches y R for unstable configuration, 
dRf * 


1 Y R x and (Y R q - M > 0). Equation (82) 


Therefore — r~ vanishes when J 
dx 

x 

suggests a law of R 1 evolution. The total translational resistance moment 
C x r 

Rj = + R^ (see (46a)) displays a behavior similar to R^ due to the con- 


stancy of R- . 
* 1 



Since y is not a constant in 
the two parameter model, the crack 
growth rate vs. (equation (76)) 
essentially differs from that of the 
single parameter model. The solid 
lines I and II in Fig. 12a represent 
the crack growth rate (£ vs. J) 
according to the single parameter 
model with two different values of 

„ , n (1 > . »< 2 > 

R^, i.e., R^ < R^ , respective- 

ly. Evolution of R^ (Fig. 12b) 
yields a transition from the curve I 
to the curve II (Fig. 12a). Various 
ways of the transition depending on 
R^ evolution are shown by dotted 
lines for monotonic transition and 
dashed line for a transition with 
crack deceleration (Fig. 12). 

Applications of the two para- 


Fig. 12 


meter model to slow crack layer 


I 
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propagation in unstable configuration of PS sheet under fatigue loading was 
recently achieved {8]. Fig. 13 taken from (8] shows a reasonable agreement 
between the theory and the experimental data for more than four order of 
magnitudes in crack growth rate. The two fatigue tests presented in Fig. 13 
were identical except of the values of o . The higher stress yields 
faster crack propagation and shorter fatigue time which was expected. Most 
important is the fact that the critical value of J^(at which avalanche-like 
crack propagation was observed) is much smaller for higher stress . The 
micrographs of the cross-section along the trailing edges of the two active 
zones shown in Fig. 13 explain this phenomenon. The lower stress produces 
more dense damage , i.e., larger translational resistance moment which 
consequently leads to larger J^. 

IX. Material Toughness Characterization in Two Parameter Model 

The necessary condition of crack layer instability (see sections VII 
and VIII) is 

J 1 “ Y * R lc (83) 

c r 

where ■ R^ + R^(t) according to (46a). Wien the sufficient condition of 

3 J 

instability (— - > 0) is met, as it happens for various loading conditions, 

die 

the condition (83) expresses the only requirement for the critical state. 
Following the conventional symbolism we introduce J^ c as a critical 
value of J^. Then using (83) and (46a) one can write 

J c “ Y *o + Y * R l (t) (84) 

* * c 

where the first term Y q « Y is a Griffith's type energy associated with 
either the crack surfaces or, in a more general sense, a near surface layer 
of intensive damage (a core of damage). During crack layer growth Y q 
assumingly remains constant. The second term describes the loading history 
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dependency of and is associated with a damage dissimination accompanying 
.the crack growth* 


Equation (84) presents J as a 



linear function of (Fig. 14). The 
slope and the intersect in Fig. 14 

* * 

gives rise to y and y 9 respective- 
ly. As shown in Figs. 15, 16, 17 
the linear relationship has 

been demonstrated by the results on 
stainless steel 11], polycarbonate 
and polystyrene, respectively [7]. 



As it was suggested [12], Y appears to be a constant quantity of the 
same order of magnitude as the latent energy of a phase transition for the 
material considered [1,7]. 

X. Conclusion 

The model described above suggests three independent parameters Y* , 

o 

Y and Rj to characterize material toughness. Two of them Y* and Y* are 
material constants, reflecting the mode of damage (microstructural fea- 
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tures) . The third is a history dependent parameter. Therefore, a complete 
toughness characterization requires establishment of the constitutive 
equation for R^. 

The model described is in a good agreement with available experimental 
data. At the same time the limitations of the model are obvious. It does 
not predict the crack trajectory (we assumed a rectilinear path), it does 
not describe the active zone shape changes, observed recently [3], etc. 
Therefore the necessity of employing the rest of the crack layer degrees of 
freedom is clear. 
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Appendix 

Historical Remarks 

J-integral was first introduced by J. Eshelby in 1951 [16] to express 
the force acting on a singularity within an elastic solid. Later, J-inte- 
gral was rederived independently by J. Sanders in 1960 [17] , G. Cherepanov 
in 1967 [18] and J. Rice in 1968 [19]. The most clear and popular inter- 
pretation of J as the energy release rate with respect to crack length has 
been done by J. Rice. W. Gunter [20] (1962) and later J. Knowles and E. 
Sternbert [21] (1972) applied Noether's theorem to elastostatics and ob- 
tained three path-independent for a linear, homogeneous and isotropic 
medium integrals J, L and M associated with translational, rotational and 
expansional invariances. These integrals express the general conservation 
laws of elastostatics. Shortly after J. Knowles and E. Sternberg' publica- 
tion, the physical interpretation of L and M integrals was discussed by B. 
Budiansky and J. Rice [22]. The same J, L and M integrals (and additional 
one N) appeared in the crack layer theory as active parts of thermodynamic 
forces reciprocal to crack layer extension rotation, expansion and distor- 
tion ([5], [6] (1978)). In recent publication of S. Aoki, K. Kishimoto and 
M. Sakata [23] a generalization of J, L and M integrals is proposed for 
cases in which plastic deformation, body forces, thermal strains may exist. 
The formalism used in [23] is very similar to that in [5,6]. Path indepen- 
dent integrals for inelastic materials are in details discussed by Stone- 
sifer and Atluri [24]. 


1. Report No. 


2. Government Accession No. 


3. Recipient's Catalog No. 


NASA CR-174634 

4. Title and Subtitle 


5. Report Date 


Crack Layer Theory 


March 1984 

6. Performing Organization Code 


7. Authors) 


8. Performing Organization Report No 


None 


A. Chudnovsky 


9. Performing Organization Name and Address 

Case Western Reserve University 
Dept, of Civil Engineering 
Cleveland, Ohio 44106 


12. Sponsoring Agency Name and Address 


10. Work Unit No. 


11. Contract or Grant No. 

NAG 3-23 

13. Type of Report and Period Covered 

Contractor Report 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. Supplementary Notes 


14. Sponsoring Agency Code 


505-33-22 


Project Manager, Bernard Gross, Structures and Mechanical Technology 
Division, NASA Lewis Research Center, Cleveland, Ohio 44135. 


16. Abstract 

Studies of crack propagation and stability have developed in two main directions: One related to 

materials science, which concerns studies of the hierarchy of microdefects, their nucleation, inter- 
action and development in association with propagation of a main crack. Another, based on continuum 
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bridge these approaches we introduce a damage parameter in addition to conventional parameters of 
continuum mechanics and consider a crack surrounded by an array of microdefects within the continuum 
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layer." Crack layer propagation is an irreversible process. Hence, the general framework of the 
thermodynamics of irreversible processes have been employed to identify the driving forces (causes) 
and to derive the constitutive equation of CL propagation, that is, the relationship between the 
rates of the crack growth and damage dissimenation from one side and the conjugated themodynamic 
forces from another. The proposed law of CL propagation is in good agreement with the experimental 
data on fatigue CL propagation in various materials. The theory also elaborates material toughness 
characterization. It proposes the following relationship between the critical energy release rate 
Jir (widely used as a material toughness parameter) and the damage dissemination characteristics 
R<j (which is called translational resistance moment): 

3, r R 1 0) 

1C o 1 1 

Here 7* is a Griff i th 1 s type energy associated with either the crack surfaces or, in a more 
general sense, a near surface layer of intensive damage (a core of damage); 7* is the specific 
enthalpy of damage. is measured on the basis of the fracture mechanics methods, 7 0 , 

7 and R] results from the materials science technique. Thus, equation (1) suggests a practi- 
cal link between the two approaches to crack stability analysis. 
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